Stringent limits on quark star masses due to the chiral transition temperature 
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Recent finite-temperature QCD lattice data are analyzed 
within a quasiparticle model, and extrapolated to nonzero 
chemical potential. Determined by the chiral transition tem- 
perature, the resulting equation of state of charge neutral, 
/3-stable deconfined matter limits the mass and the size of 
quark stars. 
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Soon after the discovery of quarks and gluons as the 
constituents of hadrons, the possible existence of quark 
cores in neutron stars has been conjectured Q. Even 
before establishing quantum chromodynamics (QCD) as 
fundamental theory of strong interaction, the notion of 
quark matter in stars had been introduced by various au- 
thors 1^ . Since that time many investigations have been 
performed to elucidate various consequences of this hy- 
pothesis. As one of the most fascinating conjectures, the 
existence of self-bound quark matter and pure quark 
stars (as surveyed in |^) was proposed. Other investi- 
gations are devoted to rapidly rotating quark stars 
the cooling behavior of neutron stars with quark cores, 
and the mass-radius relation which is now becoming ac- 
cessible to observation. If star matter exhibits a strongly 
first-order phase transition above nuclear density, a third 
stable family ||^ of ultra-dense, cold stars can appear be- 
yond the branches of white dwarfs and ordinary neutron 
stars, thus creating so-called twin stars. The necessary 
discontinuity in the equation of state (EoS) may be real- 
ized in nature at the chiral transition |Q , but also other 
phase changes of nuclear matter may allow for twin stars 
1^ . The transition of a neutron star to the denser config- 
uration, triggered e.g. by accretion, can be accompanied 
by an ejection of the outermost layers This may re- 
sult in specific supernova events [pyiC|]. 

While different scenarios for these phenomena have 
been explored, most studies suffer from uncertainties of 
the EoS of strongly interacting matter. Predictions of 
the hadronic EoS, which can be adjusted to properties 
of nuclear matter, become increasingly model dependent 
when extrapolated above the nuclear saturation density. 
For deconfined matter, which we will consider in the fol- 
lowing, the bag model EoS is often applied 0], which 
is, however, in conflict with thermodynamic QCD lattice 
data. Improved approaches take into account medium 
modifications of the quark masses or rely on pertur- 
bation theory [T^ , but also introduce a priori undeter- 
mined model parameters. 



Here we are going to describe deconfined matter within 
a quasiparticle picture [ p^Jl4| . The approach is tested 
in the nonperturbative regime by comparison to finite- 
temperature QCD lattice data, which completely fixes 
the intrinsic parameters. The model yields in a straight- 
forward way the EoS at finite chemical potential and van- 
ishing temperature, which is then employed to calculate 
bulk properties of quark stars. 

Our starting point is the thermodynamic potential in 
the form p{T,fi) — X^ift ~ -^(np, where the contribu- 
tion pi{T, fii(fi);mf) of the thermodynamically relevant 
partons is given by the pressure of an ideal Bose or Fermi 
gas of quasiparticles with mass squared mf = ttiq ^ -f 11* . 
TOo,i is the rest mass of parton species i, and 11* its 
asymptotic self-energy at the light cone (depending on 
the temperature T and the chemical potential /x). Per- 
turbatively, this ansatz is motivated by the fact that 
in the interacting system the position of the poles of 
the transverse gluon propagator and the quark particle- 
excitation is, for the dominating hard momenta of the 
order T, fi, described approximately by mass shells. The 
longitudinal gluon modes and the quark hole-excitations, 
on the other hand, are overdamped in this region of 
the phase space. The function -8(11*) is determined 
from a thermodynamical self-consistency condition, via 
dB/dn* = dpj{T,n,-m'])/dm]. 

This phenomenological approach reproduces the lead- 
ing order perturbative result for the QCD pressure Q. 
Since it partly resums terms of higher order, it can be 
supposed to be applicable at larger coupling strength. 
(This assertion was also made in Refs. where the 
pressure or the entropy have been approximated by re- 
summing the hard thermal loop self-energies.) Indeed, 
parameterizing the effective coupling by 
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finite-temperature QCD lattice data can be reproduced 
quantitatively, even in the close vicinity of the chiral tran- 
sition temperature Tc, by adjusting and A as demon- 
strated in for Nf = 2 and 4 light flavors, and in 
the quenched Hmit. Based on that, only estimates could 
be provided for the relevant case of three quark flavors 
with physical masses. Now, the predictions of the quasi- 
particle model can be considerably refined by systemat- 
ically analyzing the recent lattice calculations ||lq] with 
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Nf = 2, 3 and Nf = 2 + 1. First of all, we emphasize the 
fact that the pressure, as a function of T/Tc and scaled 
by the free value, displays a striking universal behavior 
for different numbers of flavors This observation 

supports the proposed quasiparticle description. 

In Ref . |l6) , the pressure p{T) has been calculated on 
the lattice for two and three flavors with a quark (rest) 
mass m}^^^ = 0.4 T. These data can be reproduced by 
the quasiparticle model assuming either mo,g — OAT (fit 
1 in Fig. 0a for Nf = 2) or mo,, = (fit 2); it natu- 
rally appears that the first option yields a slightly better 
fit. To estimate the effect of the temperature dependent 
quark mass m^'*'*, we also calculated the pressure within 
the quasiparticle model assuming mon = and with A 
and Ts from fit 1 (dotted line in Fig. |^a). The resulting 
variation of the order of 5% is within the uncertainty of 
the lattice data. To extrapolate the lattice data to the 
chiral limit, and to correct for finite size effects, it has 
been suggested in to scale, for 2Tc ^ T ^ 4Tc, the 
'raw' data by a factor of 1.15. Assuming, as a conserva- 
tive estimate, the validity of this extrapolation also for 
Tc ^ T ^ 2Tc (as relevant below), it can be described 
in that interval with the parameters of fit 3. For larger 
temperatures, this fit is still within the estimated errors. 
The data for Nf = 3 are analyzed correspondingly, and 
the fit parameters are summarized in Table 1. 

Encouraged by the convincing agreement of the quasi- 
particle model with the finite-temperature lattice data in 
the nonperturbative regime, it can be further extended to 
nonzero chemical potential which is not yet accessible to 
lattice calculations. This only assumes that the underly- 
ing quasiparticle structure does not change^. As shown in 
iQ, Maxwell's relation imposes a partial differential flow 
equation for the effective coupling G^{T, /x), which allows 
to calculate the pressure in parts of the T, /i plane if the 
coupling is known at finite T and /i = 0. From the pres- 
sure p, the entropy and particle densities s and n follow 
by differentiating with respect to T and /i, respectively, 
and the energy density is obtained by e = sT + nfi — p. 

The above quasiparticle fits in which mo,q do not de- 
pend on T are mapped to finite quark chemical potential 
fiq — fi and to zero temperature. The difference of 15% 
between the fits of the 'raw' lattice data and the esti- 
mated continuum extrapolation is also apparent at /i 7^ 0. 
Remarkably, however, the relation e(p) at T = turns 
out to be rather similar for both fits. In particular, in 
both cases the energy densities at vanishing pressure are 
approximately equal. For p ^ 15 T^, both curves can be 
parameterized within 2% accuracy by a linear function 
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^ For the equation of state, the phenomenon of color super- 
conductivity plays a less important role since it occurs on the 
Fermi surface. The effect is of the order (A//^)^, where the 
gap A is small compared to p 



Oh 




fit (a) 
fit (b) 



T/ 



FIG. 1. The pressure for Nf = 2 (upper panel) and 2-1-1 
(lower panel) at ^ = 0. The full symbols represent the lattice 
data Ijlql, the open symbols are estimates for the continuum 
extrapolation. The quasiparticle fits and the dotted curve are 
explained in the text. The parameters [A, —Ts/Tc] for Nf = 2 
are according to Table 1, and [3.5, 0.60] for fit (a) and [6.6, 
0.78] for fit (b), respectively. 



e{p) — AB + ap . 
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The parameter B is related to the bag constant in the 
bag model EoS, while a slope parameter a different from 
abae — 3 indicates nontrivial interaction effects. For the 
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and a[2] ~ 3.4 . . . 3.7 
the smaller one fit 
« 15 . . . 16.5 and 



3.8. The parameters are approximately 
= 2 and Nf — 3, which is, given the afore- 



3). For Nf 
Q![3] « 3.4 . . 
equal for Nf 

mentioned universality of the scaled pressure for different 
numbers of flavors, not unexpected. Nonetheless, we em- 
phasize again that in both cases the EoS at T = in the 
form (^ is to a noteworthy degree insensitive (in contrast 
to the model parameters A and Tg) to the remaining un- 
certainties of the lattice data. 

With these remarks in mind, we turn to the discussion 
of the 2+1 flavor case. The lattice calculation [M as- 
sumed to}***' — OAT for the light flavors, and to'^ = T 
for the heavy one. As shown in flt (a) in Fig. |l|b, the data 
can be described in the quasiparticle model assuming the 
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same temperature dependent quark rest masses as on the 
lattice. For an estimate of the continuum result and in 
order to extrapolate to physical quark masses, we scale 
the lattice data again by a factor of 1.15 (motivated by 
the universality of the pressure for different flavor num- 
bers and by the observed weak dependence on the quark 
rest masses). The scaled data can be reproduced by the 
quasiparticle model assuming toq,; = for the light fla- 
vors I — u,d, and mo,s — Tc for the strange flavor, see fit 
(b) in Fig |l|b. 

We now apply our model, with the parameters of fit 
(b) , for the case of neutral quark matter in /3 equilibrium 
with leptons. The contribution of the leptons is approx- 
imated by a free electron gas, so the chemical potentials 
are related by fi — ^d,s = ^J■u + Me- The electron chemi- 
cal potential Heil^-) is determined by the requirement of 
electric neutrality. The pressure, as well as the quark 
and energy densities at T = 0, are depicted in Fig. ||. 
(With respect to |l^] we mention that these results are 
not sensitive to the lepton component.) 

The pressure of the deconfined phase vanishes at a 
chemical potential Hc ~ 3.5 Tc, which provides a lower 
bound for the transition to the chirally broken phase. 
Lattice calculations determine Tc « 150 . . . 170 MeV. 
The uncertainty has, apart from the scaling, only a 
marginal effect on the thermodynamic properties shown 
in Fig. ^ Assuming Tc — 160 MeV in the following, the 
net quark density at fic — 556 MeV is approximately 
five times the quark density po — 3 ■ 0.17/fm'^ in nuclear 
matter. This corresponds to an energy per baryonic de- 
gree of freedom of about 1.6 GeV. Consequently, strange 
quark matter is not stable according to our model. On 
the other hand, the value of nq(fic) imposes a restriction 
on the EoS of hadronic matter, which still exhibits large 
theoretical uncertainties for higher densities. 
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FIG. 2. The pressure, the strange and the net quark den- 
sity, and the energy density for neutral and /3-stable decon- 
fined matter at T = 0, scaled by the Stefan-Boltzmann limits, 
where all masses can be neglected and fie vanishes. 
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FIG. 3. The relation e{p) from Fig. 2. 

The relation e(p) can again be parameterized by the 
linear function (||) (see Fig. ||) with a « 3.8 and B^/^ w 
228 MeV, which is considerably larger than the typical 
value of 145 MeV often used in the bag model. 

It is interesting to compare our results to those of 
the perturbative approach fl^ , where also an almost 
linear relation e(p) was found. For the choice of the 
renormalization scale A — 2/i, an effective bag constant 
« 210 MeV (at a = 3 fixed) and similar quark densities 
were obtained despite of a somewhat smaller value of /ic. 

The structure of nonrotating stars is determined by the 
Tolman-Oppenheimer-Volkov (TOV) equation (cf. [^), 
in which the EoS of the star matter enters only in the 
form e{p). In the following, we use the above results 
and consider stars composed entirely of neutral /3-stable 
deconfined matter. We first recall a property of the TOV 
equation: If e and p are scaled by a factor of k, the mass 
and length scales change by a factor of k^^/^. With the 
EoS (2), therefore, the mass and radius of the star are 
oc B^^/"^ 1^. Increasing values of a lead to a softening 
of the EoS, which results in a larger pressure gradient in 
the star and to more compact configurations. 

The mass-radius relation according to our EoS is shown 
in Fig. ^. Compared to results based on the conventional 
bag model EoS, which predict a maximum mass up to 
1.5Msun and a corresponding radius R « 10 km similar 
to the bulk properties of neutron stars, we find configu- 
rations which are smaller and lighter by, approximately, 
a factor of two (rotational effects may increase the max- 
imum mass up to 30% j^). This general statement is 
robust in our approach. The main effect is due to the 
large value of B, whereas the change a = 3 ^ 3.8 re- 
duces the maximum mass and radius only by about 15%. 

As pointed out in Refs. ||,^,|l2|, an EoS allowing for a 
strong first-order transition to hadronic matter, as ours, 
can lead to a second stable branch in the mass-radius 
relation nearby the neutron star peak, thus permitting 
twin stars. These features, however, are sensitive to 



3 



e 

CO 




RfkmJ 

FIG. 4. The mass-radius relation of quark stars for the EoS 
(2) with B^'''* — 228 MeV and parameters a around the fitted 
value. Left to the maximum mass, the upper branches become 
unstable. 

details of the hadronic EoS and of the phase transition 
and deserve separate investigations. 

In summary, our quasiparticle analysis of recent finite- 
temperature QCD lattice data allows to predict the 
EoS at nonzero chemical potential. At T = 0, the rela- 
tion e{p) is found to be almost linear and insensitive to 
the remaining uncertainties of the lattice data. For QCD 
with physical quark masses, the resulting energy densi- 
ties are large and exclude self-bound quark nuggets or 
strangelets. We find stringent limits on the mass and the 
size of strange quark stars: Assuming Tc — 160 MeV, the 
estimated maximum mass is Mmax ~ O.TMsun, with a 
corresponding radius of i?max < 5 km; these values scale 
as T~^. Such small but dense objects are of interest with 
respect to MACHO events |l|j2l|]. 
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fit 1 


fit 2 


fit 3 


Nf = 2 


[4.4, 0.73] 


[6.0, 0.80] 


[11.5, 0.90] 


^/=3 


[5.3, 0.73] 


[7.8, 0.80] 


[15.2, 0.88] 



TABLE L The parameters [A, ~Ts/Tc\ in the effective cou- 
pling (0), fitted to the lattice data [ ]lq ]. 
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